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I. INTRODUCTION 


Precision tests of QED in atoms were first carried out for hydrogen, and have been 
extended to other one-electron systems such as positronium and muonium, as well as helium 
lithium Pi, and even beryllium a. all l.ese caaea Ihe auclea. cKaa.e Z .a low. ao 
the basic expansion parameter of bound state QED, Za, is a small quantity. For this 
reason techniques in which the smallness of this parameter is exploited have been rehned 
in sophistication over the years, culminating in the present widespread use of effective field 
theories such as NRQED ^ and the effective Hamiltonian method However, at the 
same time experiments of increasing precision have been carried out on both highly charged 
hydrogenlike ions and also ions with more electrons, where as an example of the accuracy 
achieved at the highest Z we note the recent determination P of the 2pi/2 — 2si/2 transition 
energy in lithiumlike uranium. 


= 280.645(15) eV. 


( 1 ) 


As the expansion parameter Za is no longer small in this case, techniques in which an 
expansion in it is avoided are necessary. In the non-recoil limit, in which the nuclear mass 
is taken to inhnity. Furry representation QED P allows a systematic Feynman diagram 
based treatment of highly charged ions. A central structure in this approach is the electron 
propagator in a Coulomb held, the Dirac-Coulomb propagator, which satishes the equation 


A -r^r 7 o ^7 ■ V - m 
. \x\ 


SF{x,y;E) =6^ {x - y) 


( 2 ) 


As hrst pointed out by Wichmann and Kroll for the vacuum polarization, and by Brown, 
hanger and Schaefer jof for the self energy, treating this propagator exactly using numerical 
methods allows a determination of the Lamb shift that automatically accounts for all orders 
of an expansion m Za. 

When applied to lithiumlike uranium, use of this propagator gives a one-loop Lamb shift 
contribution (including screening corrections) to the 2pi/2 — 2si/2 splitting of -41.793 eV, 
which when combined with the nonradiative energy shift of 322.231 eV leaves a 0.207 eV 
discrepancy with experiment. This can be used to infer the two- loop Lamb shift, which has 
recently been calculated for the ground state of hydrogenic ions |l0|, but before this can be 
done recoil terms, the subject we wish to address in this paper, must be reliably calculated. 
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Recoil effects have of course been treated for low Z atoms, but again the techniques 
cannot be directly extended to high Z ions. The general level of treatment of these small 
corrections in this latter case is to scale the overall energies by a factor of /i/m, where the 
reduced mass /i is dehned in the usual way in terms of the electron mass m and the nuclear 
mass M, 


vnM , , 

^ m+M ^^ 

This has a relatively small effect for the transition we are discussing, amounting to only 

-0.006 eV, below the experimental error. A larger effect comes from the mass-polarization 

operator, 

- ■ (4) 


Hmp — X! 

i<j 


Pi ■ Pj 


M 


When this term is evaluated in a realistic potential it contributes -0.081 eV, for a total lowest 
order effect of -0.087 eV, 42 percent of the discrepancy. To accurately infer the two-loop 
Lamb shift, a more sophisticated treatment of recoil effects is clearly needed. The only 


such treatment we are aware of is that given by Shabaev [11| and collaborators |l^, who 
in fact hnd significant corrections to the above result. The present paper is intended to lay 
the groundwork for an alternative approach. We will restrict our attention to the hydrogen 
isoelectronic sequence, and in addition restrict our attention to diagrams that contribute 
in the low Z case to order m'^{ZaY/M, leaving the treatment of a more complete set of 
diagrams, along with the treatment of the many-electron problem, for a later paper. 

While considerable progress has been made in QED in recent years with the use of effec¬ 
tive held theories, these rely on expanding around the nonrelativistic Schrodinger equation, 
which, as just discussed, is not appropriate for highly charged ions. However, the Bethe- 


Salpeter formalism 


Q, i 


introduced hrst to treat the binding of the deuteron and shortly 


afterwards applied to the atomic problem by Salpeter [l^ , allows the problem to be treated 
in a systematic manner. However, this equation is famously difficult to apply, and most 
applications rely on expanding around the nonrelativistic limit, which we wish to avoid. 

The treatment given by Shabaev is fairly complicated, and we wish to provide a cross 
check by introducing as simple a formalism as possible. This can be done by slightly modi¬ 


fying a formalism introduced by Lepage j^, and it is this approach we will now describe. 

The plan of the paper is to set up in the following section a three-dimensional formal¬ 
ism equivalent in rigor to the Bethe-Salpeter equation. In the next section the one and 
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two photon exchange diagrams that contribute to order m?{ZaY/M will be evaluated in 
Coulomb gauge, and their nonrelativistic limit will be taken. This will be followed by a 
NRQED treatment, and in the conclusion we will describe how a calculation relevant to 
highly charged ions can be carried out. 


II. FORMALISM 


Q 


It has been known for quite some time m that there is an arbitrary number of bound 
state equations equivalent in rigor to the original form of the Bethe-Salpeter equation, but 
that are effectively three-dimensional. Many practical calculations have used formulations 
that also incorporate the Schrodinger equation j^. However, for the problem we are 
considering a relativistic approach is demanded. We note that a fairly detailed discussion of 
a number of notational and formal issues involved with the use of three-dimensional forms 


of the Bethe-Salpeter equation is given in Ref. [1^, to which we refer the reader interested 
in more detail. 

The spectrum of many high-Z ions has been studied, and it would be impractical to 
consider the differing spins of each nucleus. For this reason we simply treat the nucleus as 
a spinless particle of charge Z|e|. The Feynman rules for the electrodynamics of a spin-0 
particle involve the coupling iZ\e\{p + p')fM for the one-photon vertex, and 2i{ZeYg^y for 
the seagull vertex. The /x = 0 component of the one-photon vertex for a nucleus close to 
mass shell will then be dominated by the factor 2M, where M is the mass of the nucleus. 
Hyperhne effects associated with nuclear spin can be treated separately. We can also model 
the hnite size of the nucleus by replacing the nuclear charge Z with a form factor Z{q‘^) if 
desired, but this will not be done in this paper. 

We now consider the truncated two-particle Green’s function for the scattering of an 
electron and nucleus. We dehne the initial and hnal electron three-momenta as k and I 
respectively, and work in the center of mass so that the corresponding nuclear momenta 
are —k and —1. For the fourth component of momentum we choose Ei -|- k^ and Ei -|- /q 
for the electron line and E 2 — ko and E 2 — Iq for the nuclear line, where Ei and E 2 will be 
chosen close to the electron and nuclear masses, and when the total center of mass energy 
E = El + E 2 is & bound state energy a pole will be present. The formalism to be described 
below in its simplest form leads to a perturbation expansion about Ei = e and E 2 = M, 
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with e equal to the Dirac bound state energy, 


= m 1 + 


2 -, - 1/2 


= m< 1 — 


\n - (j + 1 / 2 ) + ^{j + 1 / 2)2 _ ^za)y J 

2n2 2iiS b + 1/2 4nD‘2(( “) )| . 


giving a total bound state energy E = M + e = M + m — m + ■ ■ ■ However, since this does 

not incorporate the known reduced mass dependence of the nonrelativistic binding energy, 
we will instead arrange the formalism so that Ei = S = and E 2 = M + m — fi = M. 
We incorporate these energies into the four-vectors Pi = {S,0) and P 2 = (M, 0). If we also 
dehne four vectors k = {ko, k ) and I = (/q, I ), the initial electron and nuclear momenta are 
Pi + k, P 2 — k, and the hnal momenta Pi + I, P 2 — 1. The truncated two particle Green’s 
function obeys the equation 


Gt{Pi + k,P2- k; Pi + l,P2-l)= iK{Pi + k,P2- k- Pi + l,P2-l) + 

I - 0 ^,KiPi + k,P2- k; Pi + g, P2 - q)Siq)GTiPi + q, P2 - q; Pi + /, P2 - /), ( 6 ) 

where K represents all two-particle irreducible kernels and the spin-l/2-spin-O two-particle 
propagator has the form 

7 (Pi + q) — m + ie (P 2 — qY — M2 -|- ie 
For brevity in the following we will write the above as 

1; E) = ,K{k, l-,E)+ I -0-,K{k, q; E)S(q)GT(q, k E). ( 8 ) 

The main point of all simplihcations of the Bethe-Salpeter formalism is that we can 
replace the relatively complicated two-particle propagator S with a simplihed form Sq and 
write 

Grik, /; E) = iK{k, 1; E) + J -^^k, q; E)Soiq)GTiq, /; P), (9) 

which serves to dehne K through 

K{k,k,E) = K{k,k,E)+ J ^^K{k,q;E)(S(q)-Sq(q))K{q,k,E) + 
/^/ -^l<(>=’<l’E)(.S(q)-S„(q))K(q,p-,E)(S(p)-S„(p))K(p,k,E) + ... ( 10 ) 
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Our choice for S'o(g) is 


So{q) = 


7i5{qo 


M S'Jo - ^ ■ q- fi + ie 
*7rh(go) 


M 




( 11 ) 


As mentioned above, we could have chosen another form with £ and qt replaced with e and 
m, which would lead to the Dirac equation with mass m in the M ^ oo limit. Our method 
will lead to a Dirac equation with reduced mass /i in that limit. We note that this method 
of building in the reduced mass is relatively simple, in particular requiring no rescaling of 


coupling constants. The delta function we have chosen differs from that of Ref. jl 


m 


which a delta function that puts the nucleus on mass shell is chosen. While the latter choice 
has a number of advantages when Feynman gauge is used, we use Coulomb gauge in this 
calculation, and putting the nucleus on-shell is not needed. At this point we can go to a 
completely three-dimensional formalism by choosing = Iq = 0, a choice that has no effect 
on the location of the bound state poles, and which allows us to replace the four vectors k 
and I with k and 1. In this case Eq. m takes the three-dimensional form 

1 


GT{k,l]E) = iK{kJ;E) + 


2M J (27r) 


^tK(k, f, E)S„(,--)Gt(«, t E). 


( 12 ) 


One gets to a bound state equation by creating an untruncated Green’s function G{k, 1; E) 
dehned through 


G(kX E) = - 1) + ^So(k)GT(^; E)So{r) 


that satishes 


1 


G{k,l;E) = -G^So{k){27if6%k - 1) + So{k) 


1 


^tK{k,q-,E)GiqXE). 


(27r 


(13) 


(14) 


2M 2M 

While this function differs from the Bethe-Salpeter untruncated Green’s function, it has 

□ 

poles at exactly the same total energy jl9|]. In order to obtain a solvable problem, we now 


introduce the simpler equation 


Goik, I- E) = ^S,{k){2TiY5\k -l) + SYk) 


2M J (27r) 


^iK^cik,q;E)GMk,E), (15) 


which can be written, because the kernel Kic for one Coulomb photon exchange is 

AniZa ~ 

Eic = ^ —^2M7o, 

\k — q\ 


(16) 
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as 


-I /-73 1 

(£-,„-i-k-p.)G„(kJ-,E) = ^(2^)V(k-l)-ij!Za j ( 17 ) 

where we have multiplied by SQ^(k). 

In the following section we will discuss the effect of expanding G about Go, but here 
restrict our attention to the latter function, which, except for the factor of ^ multiplying 
the delta function, is precisely the momentum space form of Eq. m with the electron mass 
replaced with the reduced mass. It has the spectral representation 

Go{UE)=j: ''’flG‘\ ( 18 ) 

n ^ 


where 'ipni^) is the solution to the Dirac equation with the usual normalization factor mul¬ 
tiplied by ^and has poles when E = M + 8n = Eq. To illustrate, the ground state (g) 
wave function has energy 


Eq^ = M + fi'^ 


( 19 ) 


where 7 = — {Zay, and the form 




2M 


{jjiZa) 


V ^f(p)x-iM 


( 20 ) 


where g{p) and /(p) can be expressed in terms of the dimensionless variable q = p/{pZa), 


9{p) 

fip) 


^sin(0(l + 7))[l + gT^'^^^^' 

Zq 


N 


g3(l + 7 ) 


- gcos( 0 (l + 7 )) 
7 


[1 + g 


2 i -( 1 + 7)/2 


Here 6 = tan and the normalization factor N is 


( 21 ) 


N = 2^+Vr(i + 7 ) 


(1 + 7) 

^r(i + 27)' 


( 22 ) 


The Xk’s are two-component eigenfunctions of J^, J^, and and are labeled by k = 
=F(j -|- 1/2) for j = i El/2 and p, the quantum number corresponding to J^- The spherical 
spinors are normalized so that x^Xi when integrated over solid angle, gives one. In the 
following we adopt the convention of working with Dirac wave functions with the usual 
normalization, which we account for by multiplying a factor 1 /( 2 M) into expressions for 
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energy shifts, which always involve two Dirac wave functions. We note that in most three- 
dimensional formalisms a lowest order potential that differs from the one-Coulomb photon 
exchange kernel must be devised to obtain a Schrodinger or Dirac equation, and it in general 
has energy dependence, which leads to derivative terms: neither complication appears in the 
present formalism. With this definition of our lowest order problem, we now turn to the 
calculation of corrections to the lowest order energy, M + S. 

III. PERTURBATION EXPANSION 

Perturbative corrections to the lowest order energy can be derived by calculating the shift 
of the pole position in G. To the order we are interested in here, the shift can be shown to 
be 

E-E„ = AE, = ^I f; E) - iK,c{k, I, EMI). ( 23 ) 

Schematically we can write 

K = Kic + Km + Kccx + Kccs + Km{S - So)Km + ... ( 24 ) 

Here IC* and IT refer to one Coulomb and one transverse photon exchange, OCX is the 
crossed ladder diagram with two Coulomb photons, CCs is the seagull diagram with two 
Coulomb photons, and the last term is the leading part of the correction induced in the kernel 
by our change of propagators as given in Eq. m- The term KicSKic is an uncrossed ladder 
diagram, denoted Kcc, and it is easy to see that the term —KicSqKic is equivalent to —iCic 
when used to evaluate AEi. The net effect, illustrated in Fig. 1, is that 

K = Kit + Kcc + Kccx + Kccs + ... ( 25 ) 

Only these four terms need be considered to obtain the corrections of order m^{ZaY/M, 
and we now turn to their evaluation. 


A. One-transverse photon exchange 


The transverse photon propagator with momentum q depends on both go and g. It 
simplihes in our formalism, which forces go = 0, and Km gives the energy shift 
AirZa rd^kdH 1 


A Em — 


■ (fc-I- / )'0(/) — g . (fc -I- / )i/j'^{k)a ■ qi/j{l )], (26) 


2M 


(27r)6 (g2)2 






f 
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FIG. 1; Contributions to energy levels of hydrogenlike ions at order m?{Za)‘^/M. Initial and final 
wavefunctions are implicit. Vertices on the top (electron) line are the usual spin-1/2 QED factors 
—(where e = |e|). Vertices on the bottom (nucleus) line are those appropriate for spin-0 QED: 
iZe{p + for the one-photon vertex and 2i{ZeYg^y for the seagull vertex. Part (a) represents 
the exchange of a transverse photon; (b) represents the exchange of two Coulomb ladder photons; 
(c) represents the Coulomb-Coulomb crossed ladder; and (d) represents the Coulomb-Coulomb 
seagull. 


where q = k — 1. If we approximate the Dirac wave functions in terms of Schrodinger wave 
functions through 


^(p) = 


this simplihes to 


AEit{NR) = 


AnZa f (Rk r dH 


AmM J (27r)3 J {2ttY 
AnZa r d^k r d?l ^ 


mM J (27r)3 J (27r)3 


' <^NR{v)XKiiv) 

\ ^4>nr{p)Xki,{p) 

4^\r{^)4^nr{1) 
^\r{^)4’nr{1) 


\k + l\^ + 2ia ■ {k X 1) {k^ — 


(27) 


q 

k'^R — {k ■ ly ia ■ {k X 1) 




+ 


2 g 2 


, (28) 


where the spherical spinors are understood. This can be Fourier transformed into coordinate 
space, leading to spin independent and spin dependent operators Hji and Hso, 


Hr 

Hso 


Zcx . ^ ^ ^ 

[dij + XiXjjpiPj 


2p?r 
Za 




L ■ a. 


After working out the expectation value of Hr one hnds 


AEr = T7^{Zaf 


^ I ^ 3 

n^{2^ + l) 


m 


+ M ^ ^ ’ 


(29) 


(30) 
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where 


< Hso >= 




(31) 


2n^i{i + l){2i + 1) 

It is of course straightforward to simply use exact wavefunctions and evaluate the integral 


numerically. The results of doing this 


mensional integration program VEGAS 


or the ground state using the adaptive multidi- 


22j | are shown in Fig. 2, where the exact result is 


compared with the NR approximation. As is also typical for the nonrecoil case, signihcant 
differences that would be poorly treated with an expansion in Za arise at high Z. We note 
that a £t can be carried out, giving 


AEir(ls) = 


m?{ZaY 

M 


-1 - 1.50(l)(Za)2 + .. 


(32) 


consistent with the known iZoif behavior. 


B. Coulomb-Coulomb ladder 


The diagram which requires the greatest care is the two-Coulomb photon ladder diagram, 
as it has a binding singularity. In addition, a new feature, present in a number of loop 
diagrams when the nucleus is treated as a scalar, is poor convergence in the integration 
over the fourth component of momentum, go) when Coulomb photons are present. Coulomb 
photon propagators, being independent of go, provide no convergence, and the remaining go 
is nominally logarithmically divergent, though that divergence vanishes by symmetry. We 
regulate this near divergence by introducing a factor A^/(gg + A^) and subsequently taking 
A to inhnity. This procedure introduces a term that will be shown to cancel when a gauge 
invariant set of graphs is considered. Including this factor into the diagram of Fig. lb gives 
the energy shift 


.(47rZa)2 r d^q r d^kdH {2M - q^f 

~ * 2M J (27r)4 J (27r)6 _ ^2 _ ^2 ^ 

_ 1 ?/i(fc)7o[(T + go)7o-7-g + m]7o^(/) 

A^ + gg |fc-g|2|^_r|2 + 


The Dirac equation can be used to carry out the k and I integrations, leaving 

AF ^ f A m-gof A^ 

2mJ {27iY {M - qoY - + i6 + qi 

^(g)(^7o-7-g-h)[(^ + go)7o - 7-g + m](T7o-7-g-/r)V>(g) 

((£: + go) 2_^2 + , 5 ) 


(33) 


(34) 
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FIG. 2: Plot of /S.Et for the ground state as a function of Za in units of vr?(Za')‘^jM. The data 
points represent the result of a numerical evaluation of (HH) without approximation. The solid 
curve is the function —1 — (3/2 )(Zq;)^, which includes the first relativistic correction. We note that 
the expansion in Za works well for small values of Z but fails badly for large Z. 

where = \/q‘^ + and uiq = \/q^ + We note the ‘mismatch’ in the numerator 
between terms with qt and with m: the former come from the formalism, and the latter from 
the electron propagator in the diagram, which is not altered by the choice of formalism. We 
now carry out the go integration by closing above with Cauchy’s theorem: once this is done 
we are free to introduce a new four-vector q = {S, q). Three terms result, with the simplest 
arising from the regulator, 

^Ecci = / 7^^(<f)(^-h)7o(^-h)V’(<f)- (35) 

AM J (27r)'^ 
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This term contributes in order m?{ZaY/M, but as it will be shown to cancel we do not 
evaluate it explicitly. The other two terms are both non-recoil, with the most sensitive being 

r ^ 


^EcC2 = 


« h / ^ n9 1 

+ M -ng)jo-j-q + m]{i^-fi)Y{q) 
{s + M- 


If we use the fact that 


(M + VtqY ^ ^ 0{1/MY 

{2M){2Qq) ^ ’ 


this simplihes to 


^EcC2 = 


<fq ii{q){4- Ii)[{e + M - - 7- ^)j/i(g) 

(27)3 


(E + M - a,f - 17| 


( 36 ) 


( 37 ) 


( 38 ) 


We proceed by rearranging the interior numerator in the above as follows: 


7 °(T + M - rig) - Y ■ q + m = (^+p)+ 7°(M - fig) -h (m - p) 


'+"> + ■ '’“h} 


— (^+h) + ^|(^+h) + 7°(^ ~ 






(39) 


2 m i 

where in the last manipulation we have replaced m with p, with the difference being higher 
order in 1/M. If we now define k = q"^ — and restore the factors ^ — /i on the left and 
right of the interior numerator we get 

r2 


K 


^ m\ . . . m . . . 

I + jj) (li - 


7°(m — T) -|- 7 ■ q — 7 ' 


ol_ 

2 m 


/i). 


( 40 ) 


We further make the expansion of the denominator 


{S + M- Qgf -q^ - + 2T(M - fig) - (m^ - + 0{1/MY 

= K + ^(2m^ — qY — 2m^/M + 0(1/M^) 

= K -I- ^ \S{S‘^ — q^ — nY) — S^ + 3SnY — 2m^| -!-■■■ 
= K ^Sk — (m — SY{2m -|- T) ■ -j 


= K 


M 
1 + 


S {m — SY{2'm + S) 


M 


Mk 


+ 


(41) 
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Combining these two forms then gives 


^EcC2 — 




(27r 

m 


M 


Mk 


Mk 


+ 7-g-7°^j (^-/^)}V'(<f) • (42) 


If we label the terms in the curly brackets in (I42|l parts 1-6, the various contributions are 


A£'cc2i =< C >d 
A-Ec'C'22 


n4/ -1 1 

M ^ ^ 12n4 j ’ 


m 


>7A1 


AEcc 23 = —miZay , . , , 


m 




AEcc 24 = jf 'iniZa) j , 


-1 


TTl ( ^ 

/^Ecc2i = JJ 


Vol ™ , 


( 43 ) 

( 44 ) 

( 45 ) 

( 46 ) 

( 47 ) 


A-Bcc26 /^f’MZa) + 1) 


(48) 


To cancel the factor of k in the denominator, we note that its nonrelativistic limit is —q^ 
{qiZaY/n^, so that the Schrodinger equation reads 

d?p 1 


-K 


2 , -4>NR{q) = dvrZa . ^ 

2/i J (27r)-^ |<? — P I 


(1>nr{p)- 


(49) 


The strategy for the last three terms is then to ‘undo’ the Dirac equation so that the Coulomb 
potential is explicitly present, take the nonrelativistic limit as in Eq. (EZD, and then use the 
Schrodinger equation. We illustrate this with AEcc 25 , 
f cPq m 


AE, 


CC25 


(27r)6 J |^_ 1 12 2mM 
The total contribution of term CC2 is then 

AEcC2 =< V >Dirac ^ ^ 


V’(?)(#-P)7 W(#-P)V’(?) 

^(^)7o7 W7oV’(0 

(/>^nr(^)(^- l + ^-ka- q)(t)NR{l) 
0^7Vi?(?) [q-l + iB-{qxl )) 0jvi?(^ )• (50) 


J (27r)^ Mr 

{AnZay f (Ekd^qdH m 

(27r)9 J |fc-g|2|g-r|2MK 

(dTrZa)^ f d^kd?qd?l m 


(27r)9 J \k-^\2\^-f\2 2MK 

dvrZa; f d?qd?l 1 


n3(2£+l)J M 


m 


—yr < 2i^so > , 


(51) 
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where the formalism subtracts off the first, non-recoil term. It is of interest that had we used 
the formalism with S ^ e, M ^ M mentioned above, the cancellation, while still removing 
the nonrecoil term, would leave a contribution that can be shown to start with the term 
^ 2 Mn^ , the standard reduced mass contribution to the nonrelativistic energy. As we have 
chosen to build this into our lowest order solution, the cancellation is hner, and leaves terms 
starting in order im?{ZaY/M. 

The remaining part of the CC calculation involves closing around a negative energy 
electron pole, which while leading to higher powers of Za than when the nuclear pole is 
encircled, is also nonrecoil. Its full contribution is 


^Ecc3 — 


d?q 


2M J {2ttY2uj, 


■{2M -\- S -\- UJq) 


2V’(g)(^-h) [-^g7o 




[S M UJqY 




(52) 

but we can approximate S = m, Uq = m, and flq = M in the nuclear propagator, leading to 
the simpler expression 


AEcc 3 = ^ [ - j2M + 2mf _ ( 53 ^ 

2MJ{2Tif2uq{M + 2mf-M^^^^^ ^ ^ 

We will show below that although this term is nonrecoil, starting in order m{ZaY, the 
nonrecoil part cancels with a contribution from the crossed Coulomb ladder. 


C. Crossed Coulomb diagram 


The crossed ladder (CCX) diagram of Fig. Ic is given by 


AEccx = i 


,{4:7rZaY f d^q f d?kd^l 


(2M -|- go)^ 


2M J i2TTYJ {2 ttY ^M + qoY-\q-k-l\^-M^ + zS 
1 ^(fc)7o[(T-hgo)7o-7 W + "i]7o^(0 


(54) 


|fc-gp|g-rp {{S + qoY - + iS) 

Taking a pole of the regulator term gives the same result as with the ladder, thus doubling 
AEcci- While the pole from the nuclear line enters in the entirely negligible order tjv'/M^, 
the electron pole contributes at the level of non-recoil hne structure, and is 


AE, 


CCX — 


{AnZaf f d^q r d^kdH {2M - S - Uqf 


2M J{2TifJ {2Tif [M-S-ujqY-\q-k-l\^-M^ 
1 1 


\k — q\^\q — I p 


Y{k )7o[-^^q7o -l-q + mYoYil- )• 


(55) 
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It is again legitimate to make the approximations S = m, Uq = m, and M‘^+\q—k—l p = 
and further use of the Dirac equation gives the approximation 

= g / (2^)3^ (56) 

This term can now be combined with the CCS contribution, and the nonrecoil term can 
easily be seen to cancel. There remains a nonvanishing recoil term of order m‘^{Za)^/M, 

AEccs + ^Eccx = ^ [ 20 Iso ^(g)(^-h)[-^g7o - 7-<f+^](^-h)^(g), (57) 

2M J [27r)'^2ujq 

which we keep, although beyond the order of interest we are considering here, as it has an 
interesting connection with the seagull diagram. 


D. Seagull diagram 


A novel feature of our formalism is the presence of so-called seagull graphs. In Coulomb 
gauge the seagull graph consists of a Coulomb-Coulomb (CC) term and a transverse- 
transverse term (TT), with the latter beyond our present order of interest. Again regu¬ 
larizing the go integration, the CC seagull graph (see Fig. Id) contributes 


AEccs = —2i 


.{AnZaY f d'^q f d^kdH h? 


J (27r)6 A2 gg 

_ Yjk )7o[(^ + Qoho - 7 • g + mYoYjl ) 

|fc-g|2|g-r|2 ((T+ go)2-g2-m2+ i(5) 


2M 

1 


(58) 


If we again close above to carry out the go integration the regulator term contributes 
—2AEcci- as AEqci was doubled from the crossed Coulomb diagram, this completes the 
cancellation of contributions arising from the regulator term. The other pole picks up a 
negative energy electron contribution, and gives, using the Dirac equation, 

AEccs = -2^ f ^(g)(^-h)[-^97o-7-g + ^](#-h)^(g)- (59) 

2M J {2nY2ujq 

The net result is that the only role of the seagull diagram to order m‘^{ZaY/M is in canceling 
the regulator terms from the ladder and crossed ladder, and in addition it combines with 
nY{ZaY/M terms coming from the negative energy pole terms in those diagrams. 
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IV. TOTAL AT ORDER {Zaf 


The combination of the CC, CCX, CCs, T, and -C graphs gives 

-1 


AE = 


^ ^4 / -1 1 

M ^ ^ ISn^J 


( 60 ) 


In order to find the total recoil contribution at this order, one must combine this with the 
recoil contribution from /i(/(u, j) — 1 ), which is —(m^/M)(/(n, j) — 1 ) where f{n,j) = e/m 
is defined through Eq. ©• The total recoil contribution through terms of order (Za)'^ is 


AE, 


m 


recoil — 

“ M I 2n2 
This is the known Barker and Glover result 




8 n^ 

4 /-I 


+ 


\2n^ 2n3(j + 1/2) 

for the recoil contribution at this order. 


(61) 


V. SALPETER CORRECTION 


One of the first accomplishments of a fully relativistic treatment of the two-body bound 
state problem was Salpeter’s discovery Q| that corrections of order am/M times fine struc¬ 
ture were present in hydrogenlike atoms: use of the older formalism of the Breit equation had 
found no such contributions While we are not calculating all such terms here, we show 
how they arise from one and higher Coulomb exchanges that are crossed by a transverse 
photon, illustrating with the graph of Fig. 3a. This gives rise to the somewhat complicated 
expression 


AE = 


(AnZa)^ 

r d^k 

r dH 

r d?p d?p' 

(2M -1- 2/co — Iq) 

2M J 

' (27r)4j 

' (27r)W 

' (27r)6 

(M + fco)2 -\k-p\^ - +id 


{2M ko — Iq) 


(^-2p),(5p-^) 1 


[(M + ko- hf - - /12 - M2 + i5] 


kn — k'^ + id 


:7o- 


\l — p |2 p'|2 

—7oV’(p)(62) 


(E -F /co)7o - ^ ■ {k + p') - m + id {£ + /o)7o - ■ I - m + id 

We wish to show that this expression is dominated by a term that contains part of the 
Dirac-Coulomb propagator of Eq. ©. Specihcally, if we expand the momentum space form 
of that equation in powers of the Coulomb potential, the term involving one potential is 


Sp {P,P'',E) = 


—AnZa 


7o- 


E'yo — ^ ■ p' — m \p' — p |2 Ejo — 7 ■ p — m’ 


(63) 
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vz? 


(a) 


(b) 


FIG. 3: Graphs contributing to the Salpeter correction of order m?{Za)^/M. In graph (a) a 
transverse photon crosses two Coulomb photons. In graph (b) the transverse photon crosses a 
single Coulomb photon, with a Coulomb ladder photon on the side. There are two graphs like (b) 
since the ladder photon can be on either side. 


where the distinction between m and /i is dropped in this section as the graph being con¬ 
sidered has a factor m/M. To see how this arises from the graph of Fig. 3a, we hrst close 
the ko contour above and take the transverse photon pole to get 


AE = i 


. (dvrZa)^ r d?k 1 


AM 



f d^p d^p' ( 2 M — 2k — 1 

0 ) 

1 (27r)4j 

1 

1 

1 

CM 

P - M 2 


( 2 M -k-lo) 


[{M-k 


lof- -M2 + i5] 

1 


:-2p), 




1 


\l — p |2 \k — I A- p'P 

(64) 


--7o-- 

{S - k)'yo - A ■ {k + p') - m {S + lo)'lo - j ■ I - m + i6 

Because we are dropping terms of order /M"^, the hrst nuclear denominator simplihes to 

—2Mk, and further carrying out the Iq integration by closing above gives a term from the 

second nuclear denominator that forces Zq = —k, giving 


AE = 


{AnZa)^ f d^k 1 f dH f d^pd^p' 


AM 

^{p)%- 


(27r)^ J (2tiY 
1 


(27r)6 

-7o- 


-2p)j 




) \l — p p \k — I + p'P 
■7o^(p)- (65) 


{8 — k)'yo — 7 ■ (Zc + p') — m {8 — k)'yo — 7 ■ / — m 
The Salpeter correction is associated with the region of integration in which k ~ m{ZaY 
rather than the normal m{Za). In this region one can approximate k + p' = p' and 
k — I + p' = —I + p' in the above. Using Eq. (P|) then allows us to write 


AE = 


(ATrZa)^ r d?k 1 r dH r d?pd^p' 

AM J (27r)^ J (27r)^ J (27r)® 

^{p'hiS]f^{p',k,8 - /c)7oV^(p). 


-2p), 


5ij - 


kikj' 


kP‘ ) \l — p\ 


( 66 ) 
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The same kinds of argument apply for any number of Coulomb exchanges, allowing the 
replacement of with Sp- Care is required, however, for the hrst term of the expansion of 
Sp, which requires an ultraviolet cutoff in the k integration because of the approximations 
we have made. However, this term is hnite when simply treated as a one loop diagram. The 
actual calculation of the Salpeter correction would involve evaluating the one loop diagram 
without approximation and then evaluating the above expression and higher Coulomb ex¬ 
changes by using Sp —So, a technique that is standard in self-energy calculations. Replacing 
in the above with the spectral representation of Sp then gives 


AE = 


{AirZay r (Ek 1 r d^l r d^pd^p' 

AM J (27r)3P y (27r)3 J (27r)6 ' 

)7O'0(P ) 


■2p)j I 


kikj 


k? ) \l — p\ 


E - Em - k 

Replacing —2pj with —2(p — l)j leads to the integral 

f dSp d^l 


(67) 


J ( 27 r )8 

/ X ' — 

d^X^tljm{x)'-io'^{.x) 

= {E - Em) < m\pj\Q > 


( 68 ) 


where E is the Dirac energy of the state V’ of interest, here taken to be the ground state. 
Eqs. (iniD and (jEHl) lead to the relativistic generalization of Salpeter’s Q Eq. (45). The 
replacement of —2pj by —2(p — l)j arises from consideration of the reducible graph shown 
in Fig. 3b, which comes from the formalism when three-photon exchange is considered, as 
described in Ref. j^. In this diagram there are two nuclear propagators that depend on 
lo, with one of them leading to the replacement mentioned above, and the other to a bound 
state singularity canceled by the formalism. 


VI. NRQED CALCULATION OF THE ENERGY SHIFT TO ORDER (Za)^ IN¬ 
CLUDING RECOIL 


We base our expression of NRQED on the work of Kinoshita and Nio 
Lagrangian for a particle of spin-1/2 and one of spin-0 has the form 

^ , J ^ 3^ 3^ ea-B 

^ - ^ r * + ^ + 


2J. The NRQED 
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( 69 ) 


e(D-E-E-D) iea ■ (D X E - E X D) 
+CD- -^ + Cs - 


f 


Sm? 


>il) 


+ C 


EM 


where Dt = dt + ieA^ and D = d — ieA for the electron, and similarly but with e ^ —Ze for 
the nucleus. Electromagnetism is described by the usual Lagrangian Cem = (— 
and we use Coulomb gauge in our calculations. In terms of Feynman rules, we build on those 
given by Kinoshita and Nio in their Fig. 3. The new rules include a Coulomb vertex for the 
nucleus —Ze, a dipole vertex for the nucleus Ze{p' +p)/(2M), and a propagator for the 
nucleus {E — p^/(2M) + ie)~^. The rule is to multiply all propagators by i, all vertices by 
—i, and to include an overall factor of i when calculating energy shifts. Loop integrals are 
done over all momenta with measure d'^k/{2E)^. The Bethe-Salpeter equation for NRQED 
(with lowest order propagators and vertices) is exactly the Schrodinger-Coulomb equation 
with reduced mass, which is also described in Kinoshita and Nio. We use the symbol '^{p) 
for the NRQED Bethe-Salpeter wave function. For example, the ground state wave function 
is 


T(p) = (27r)5(po)V’(p), 

(70) 

where 


(^2 + ^2)2^-D 

(71) 

with ^ = pZa and, for example, 0)- 



The relevant graphs for the calculation of energies up to 0{{Za)‘^) and including recoil 
(to hrst order: m/M) are shown in Fig. 0] In the calculation of bound state NRQED graphs 
we take the electron line to enter with momentum (To,p) and the nucleus with momentum 
(0, —p) where To = — p{Z/{2n^) is the Bohr energy level. 

The crossed Coulomb ladder (Fig. 4a) has the form 


Am ■ f dH (Pp , + , . . 

= ' 1 


^0 + To — — ie 


'-ie)i){p) 


X- 


{q-lf {l-pf 
The poles of the Iq integral 

f dlo 1 


(iZe) 


j {i-p-qV I .• 

‘0 - m— + * 


2m 

-{iZe). 


(72) 


le 




= 0 . 


(73) 


2m 


le 
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FIG. 4: NRQED contributions to energies at order m?{Za)^/M. Graph (a) is the two-Goulomb 
crossed ladder, (b) represents the relativistic kinetic energy correction to the electron line, (c) is 
the spin-orbit correction to the electron line (with Goulomb photon exchange), (d) is the Darwin 
correction to the electron line (with Goulomb photon exchange), (e) shows transverse photon 
exchange with dipole vertices on both electron and nuclear lines, and (f) represents transverse 
photon exchange with a Fermi vertex on the electron line and a dipole vertex on the nuclear line. 


are both on the same side of the real axis. It follows that the Iq integral vanishes, as does 
the crossed Conlomb ladder contribntion: AEccx = 0. 

The relativistic kinetic energy correction (Fig. 4b) is 





(74) 


where 


(Hk) = 


/'-Ph ! 3 




i’ip) = p{Za) 


{ 27 r)^^ J ' \ 8 n^ { 2 i + l)n^^ 

The spin-orbit correction to the electron line (Fig. 4c) is 

AEso = i J (? X P) ■ 


d^q d^p 

Am? J ( 27 r)® 


{q X p- a)^p{p)Vc(k) 


= (-) {Hso) 

\m/ 

- S) <«-> 


(76) 


(76) 
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where k = q — p and 


as before. 




The Darwin correction to the electron line (Fig. 4d) is 
^ ^ , fd^qd^p ie 


AEn = i 


mj 8/i^ 

2m \ , 


where 


'-M <"«>■ 


w = i^iw)t=M^ar|| 


The dipole-dipole transverse photon exchange contribution (Fig. 4e) is 


AEnn — i 


d^q d^p 
(27r)6 




—47rZQ; fd^qd^p — (p ■ ? ^ 

~ mM J (27r)6 ^ ^ 

777 

= ^ (^k) , (80) 

just as in the relativistic calculation of one transverse photon exchange. We note that the 
expectation value (Hfi) can be written as 

= 3 (Hk) + 2 (Ho) - ■ (81) 

Finally, the Fermi correction (Fig. 4f) is 

^ -47rZa r d^qd^p -ieijk{q-p)j(Tk5fn{k){p + q)n ^,^. 

^ AmM J (27r)6 ^ ^2 


2i f d?q d^p 

~ AmM J {2nf 
m / Za ^ ^ 

M \ 2m‘^r^ ^ 

The sum of all contributions is 


^^( 7)7 X p- d'iy(p)Vc(k) 


AE — AEccx + AEx + AEqo + AEn A- AEnn A- AEi 
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( 83 ) 


3m 


M 


2m 


= 1 - ITT {-»/>■)+ 1 - ^ (Hso) +[1-^ (Ho) 


M 


2m 


M 


+^{HH,)+2(Ho)--MZar) + !^^(2Hso) 

m? {ZaY 


= {Hk + Hso + Hd) 


M 8n^ 


which again is the known Barker and Glover result for the hne structure with recoil correc¬ 
tion. 


VII. CONCLUSIONS 


A form of the Bethe-Salpeter equation of particular simplicity has been introduced that 
can be applied to the entire hydrogen isoelectronic sequence. We have shown that the 
power series expansion of the IT kernel is nonperturbative at high Z, demonstrating the 
need for a complete numerical calculation for all kernels. Such a calculation has been done 
using a Green’s function formalism by Shabaev and collaborators ^|, but it is always 
desirable in QED to have checks on these complex calculations. We are presently calculating 
the remaining one-loop diagrams that enter in order m?/M{ZaY. As an indication of the 
numerical importance of these calculations for the transition discussed in the introduction, 


we note that Ref. 


12j | Ends a correction of -0.04 eV for the 2pi/2 —2s transition in hydrogenic 


uranium, to be compared to the 0.207 eV discrepancy presumably dominated by the two- 
loop Lamb shift. However, this is only part of the effect of recoil for lithiumlike uranium, 
and the question of relativistic corrections to mass polarization cannot be addressed in 
our formalism, which is strictly a two-body approach. We are presently investigating the 
relatively unexplored problem of forming many-particle generalizations of the Bethe-Salpeter 
equation that have the three-dimensional and relativistic aspects of the equation described 
in the present work. 
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